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Abstract. We prove ergodicity of the finite dimensional approximations of the three dimen- 
sional Navier-Stokes equations, driven by a random force. The forcing noise acts only on a few 
modes and some algebraic conditions on the forced modes are found that imply the ergodicity. 
The convergence rate to the unique invariant measure is shown to be exponential. 
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1. Introduction 



The uniqueness of statistical steady states for the Navier-Stokes equations is a less known 
but nevertheless important problem in the mathematical theory of turbulence. The question 
is completely open in dimension three, mainly because, due to the lack of uniqueness of the 
equations, there is no way yet to give meaning to the mathematical objects involved in the 
O ' subject. 

In the present paper the property of ergodicity is proved for the finite dimensional approxi- 
\ mations of the three-dimensional Navier-Stokes equations, driven by a random force. The same 
O | problem has been solved in two dimensions by E and Mattingly [0. 

Such result can have a qualitative interest for the statistical behaviour of an incompressible 
fluid. Indeed, if the Kolmogorov theory of turbulence is taken into account, one can believe that 
the cascade of energy, responsible of the transport of the energy through the scales, is effective 
in the inertial range so that at smaller scales only the dissipation ends up to be relevant. Hence 
the long-time statistical properties of the fluid can be sufficiently depicted by the low modes 
of the velocity field. In some sense, if the ultraviolet cut-off is sufficiently large, in order to 
capture all the important modes, the corresponding invariant measure gives the real behaviour 
of the fluid. In view of these considerations, the conclusions of the paper can give also both a 
hint and a possible starting point for the analysis of the infinite dimensional case. 

We consider a finite dimensional truncation of the three dimensional Navier-Stokes equations, 
driven by a random force, with periodic boundary conditions. The proof of ergodicity is classical 
and it is developed in two steps. Firstly we prove that the transition probability densities are 
regular, by checking that the diffusion operator is hypoelliptic (the Hormander condition). Then 
we show that the Markov process is irreducible, in the sense that each open set is visited with 
positive probability at each time. For this aim we study the associated control problem (see 
Section |6|). Irreducibility for the infinite dimensional equations was firstly proved by Flandoli 
H, under the assumption that the noise acts on all modes. 
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Both these properties, strong Feller and irreducibility, are implied by an algebraic condition 
on the set of indices corresponding to the modes forced by the noise. The condition essentially 
means that it is possible to obtain any index as a sum of some of the forced indices. One can 
see this mechanism as a geometrical realisation of the cascade of energy, since the non-linear 
term transmits the random forcing from the few forced modes to all the other modes. As an 
example we show that the algebraic condition is satisfied if the three lowest modes are forced. 

Recently, many authors have applied the techniques we have used here, such us the hypoellip- 
ticity for degenerate diffusions, or the general theory for Markov chains and Markov processes 
collected and developed by Meyn and Tweedie (one can see their book |]IU|| ). Among many oth- 
ers we quote the papers by E and Mattingly Eckmann and Hairer M, Hairer |J, Rey-Bellet 
and Thomas [1131, and some of the references therein. 



The paper is organised as follows. In the first section the main definitions are given, together 
with the statements of the main results and an outline of their proofs. The technical compu- 
tations and the precise statement of some hypotheses are postponed in the following sections. 
The aim is to give a light presentation of the main ideas, without all the technicalities, which 
are then reserved to the interested readers. 

Acknowledgements. The author wish to thanks R. Bianchini and S. Dolfi for the helpful 
bibliographical suggestions on the control theory part in Section |6| and on the algebraic part in 
Section |5], and F. Flandoli for the many helpful conversations. 

This paper is dedicated to the memory of my father, who died whilst I was writing it. 

2. The main theorem 
We consider the stochastic Navier-Stokes equations with additive noise 

du = (uAu - (u ■ V)w - VP) dt + dB t 
divw = 0, 

in the domain [0, 27r] 3 , with periodic boundary conditions, where u is the velocity field and P 
is the pressure field, and B t is a Brownian motion. As usual, the equations are projected on 
the space of divergence-free vector fields, in order to cause the pressure to disappear from the 
equations. If we write the equations in the Fourier components, we obtain the following infinite 
system of stochastic differential equations 

du k = [-u\k\ 2 u k - i ^ (k-u h ) (u y - Tj^ k )] dt + q k d$, k G Z 3 

h+l=k ' ' 

with the constraint u k ■ k = (it comes from the divergence-free condition). We have made 
some simplifying assumptions on the noise: we assume that the noise takes values in the space 
of divergence-free vector fields and that the covariance is diagonal in the Fourier components 
(the assumptions will be stated more clearly in Section ^T|) . 

In order to state the problem of the finite dimensional approximation, fix a threshold N 
and consider the finite subset of indices K-n = { k G Z 3 | |k| < N, k ^ (0,0,0) }. The finite 
dimensional system obtained is the following 

(2.1) du k = [-^|k| 2 u k - i ^2 (k ■ ^h) (^i - kji^' k)] dt + g k d$, k G fC N 

h+l=k 
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(a formal derivation is given in Section |3|). We will use the real variables r^, Sk G R 3 , where 
«k = Tk + isk, rather than the complex variables u^, so that the equations are briefly written 

as 

dri = F r i(r, s) dt + q r k dfi, . _ 

and ?k = 9k + i?k- Since n_k = Wk, the set of indices /Cat is redundant, hence we take a smaller 
set /C, which takes into account the symmetries. 

The solution (r(t),s(t)) of the above stochastic equations is a Markov process on the state 
space 

C/ = 0(i?k©5 k ), 
kefc 

where and ,5k enclose the divergence-free condition r k • k = s k • k = (see also ( |4. 1| ) and 
the following formulas). We denote by Pt the transition semigroup 

P t <p(r ,8 ) = E (r0tSo) [(p(r(t),s(t))] 

with generator 

(2-2) £ = F Q + \ (XL/ + X } 



2 
ke/c 

i=l,2,3 



a 2n 
k,i J 



where 



(2.3) F o = Y.J2 F rl]fri +F siQ-r> X k,i = ^9k^^J, X k,i = X)3k,<i^7' 
ke/c i=1 k k 3 =1 k i =1 k 

and by Pt{(r, s), •) the transition probability. 

The main assumption we take on the noise is that it acts on a small set of modes. We 
consider the set M of indices whose corresponding Fourier components are forced by the noise. 
We assume that AT is a determining set of indices, as defined in Section |, which essentially 
means that each index in /Cat can be obtained as the sum of elements of N '. In other words 
M should be an algebraic system of generators of Z 3 . In Section |5] we will give some heuristic 
justifications to such claim. As a working example, Proposition E72] shows that any set M 



containing the three indices (1, 0, 0), (0, 1, 0) and (0, 0, 1) is a determining set of indices. 

Here we are interested in stating the main result of the paper, namely the ergodicity of the 
finite dimensional approximation fl2.1|) 



Theorem 2.1. Assume that the Brownian motion B t satisfies the assumptions in Section 3T7 
and that the set M defined above is a determining set of indices. Then the system Q2.1p admits 
a unique invariant measure. 

Moreover, the unique invariant measure is supported on the whole state space or, in other 
words, it gives positive mass to each open set. 



Proof. First, we prove the existence of the invariant measure. The method is classical and based 
on the Krylov-Bogoliubov method (see for example Theorem 3.1.1 of Da Prato and Zabczyk 
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pj). The compactness follows by the following argument. Let ||"u|| 2 = YlkeK, I n k | 2 , then by Ito 
formula (using also the first property of Lemma pM|) , 



d\\u{t)f = ^(2^-F±(u)+Tr{&'fr))dt + 2^v£-q k d0 i 
= ~2u^2 |k| 2 |u k | 2 d* + 2^Hk"- qkdffi + a 2 dt, 



ke/c ke/c 
where a 2 is the variance of the Brownian motion Bu so that 



E\\u(t)\\ 2 + 2v [ \\u(s)\\ 2 ds < E\\u(0)\\ + o 2 t 
Jo 



— 2 

and by Gronwall lemma E||w(t)|| < E||w(0)|| + |-. 

Uniqueness of the invariant measure is proved by means of the Doob uniqueness theorem 
(see for example Theorem 4.2.1 of Da Prato and Zabczyk Jl|). We just need to show that the 
transition semigroup generated by the dynamics ( |2.1| ) is strongly Feller and irreducible. 

A Markov semigroup P t is strongly Feller if P t ip is bounded continuous in time and space 
when if is bounded measurable. By a theorem of Stroock WM, the transition semigroup is 



strongly Feller if the Hormander condition holds: the Lie algebra generated by the vector fields 
in (|2.3|), evaluated at each point, is the state space U. Since N is a determining set of indices, 



from Lemma ^]2| it follows that the constant vector fields of the generated Lie algebra span U. 

A Markov semigroup is irreducible if it gives positive mass to any open set for each initial 
condition and each time. It is well known (see Stroock and Varadhan ]T5| ) that irreducibility 
is true if the control problem (see equations ( |3.1D ) associated to problem (|2.1| ) is controllable. 
The last statement follows from Theorem |6.5| . 

Finally, the irreducibility property implies also that the support of the invariant measure is 
the whole state space. □ 



The next theorem shows, by means of general techniques developed in Meyn and Tweedie [ nfl , 



T2| , that the finite approximation of Navier-Stokes equations has good dissipation properties, 
strong enough to ensure the exponential mixing of the dynamics given by the Markov process. 
In order to state the result, define, for any measurable function / > 1 and any signed measure 
fi on the Borel sets of U, 



\u,\\f = sup 

!<?!</ 



g(x)n(dx) 



and set 

^M = E E « 2 + 4 2 ), (r,s)eU. 

Theorem 2.2. Under the assumptions of the previous theorem, let n be the unique invariant 
measure. Then there are positive constants C and p such that for each initial condition (r , s ) G 
U, 

2 

||P t ((r , s ), ■) - z\\ f < Ce- pt (l + V(r , s ) + t > 0, 

where f = 1 + V . 
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3. The Navier-Stokes equations in the Fourier coordinates 

In this section we derive the equations of the finite dimensional approximations of the sto- 
chastic Navier-Stokes equations, with additive noise, 

du = (isAu - (u ■ V)it - VP) dt + dB t 
divw = 0, 

in the domain [0, 2tc] 3 , with periodic boundary conditions, in the Fourier components. 

Consider the Fourier basis (e lk:E ) kgZ 3 of L 2 ([0,2n} 3 ). First, assume that the applied random 
force has zero average, so that the centre of mass of the fluid moves with constant velocity. 
Hence, without loss of generality, we can assume that 

u = P = 0. 

The projection onto the space of divergence-free vector fields is defined as 
V(ae^) = (a- ^ ■ a)e ik - = (a - ^k)e ik -, 

where | • | is the Euclidean norm in R 3 . Notice that 

div u = means k • = for each k. 

3.1. Assumptions on the noise. For the sake of simplicity, some simplifying assumptions 
will be done. First we assume that the covariance Q of the noise is diagonal in the Fourier 
basis, so that we can write 

Qv= ^(g k -^)e lk -*. 

keza 

Moreover we assume that q£ ■ k = for each k, this implies that the Brownian motion takes 
values in the space of divergence-free vector fields. The Brownian motion has finite variance 
that we denote by a 2 . We assume also that for each k, the real and the imaginary parts of 
the 3x3 matrix q^, if not zero, have rank 2. This is an assumption in the small of non- 
degeneracy, since we ask that, if a mode is forced, it is fully forced in its 4 components. As a 
first consequence of our assumptions, the operators Q and V commute. 

The main assumption of the paper is that the noise acts only on a few components, namely 
most of the matrices are zero. We define the set N C Z 3 of stochastically forced indices, 
that is the set of ks such that q^ ^ 0. 

3.2. The equation in the Fourier modes. We write 

u (t,x) = Mt)e ik ' x 
kez 3 

and, by means of the operator V, we project the equations in the space of divergence-free 
vector fields, so that the pressure disappears. We obtain the following infinite system of sto- 
chastic differential equations (see also Gallavotti 0, Chapter 2, where the author gives also a 
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interpretation of the physics of the fluid in terms of the Fourier coordinates) 

c ■ u 



du k = [-u\k\ 2 u k - i ^ (k-u h ) (it! - ^rri^k)] dt + g k 



h,iez 3 

h+l=k 

« k • k = 

where (0£)t>o are independent three-dimensional Brownian motions, and the nonlinear term 
has been obtained in the following way: 

V(u ■ V)u = iP ^ ^ (1 • M h )«ie ik ;E 
kez 3 h+i=k 



= ^ E (i-^h)h-^-^k) 



kez 3 h+i=k 1 1 

kez 3 h+i=k 1 1 

3.3. The finite dimensional approximation. Let JVeN and set 

K. N = { k £ Z 3 | k + (0, 0, 0), IkU < iV } . 

where | • \^ is the sup-norm in R 3 . We project the equation in the space spanned by (e lk-x )k 6 £jv> 
with coefficients in R 3 , and for this aim we set 

u(t,x) = ^2 u k e±kx - 
ke/Cjv 

The equation in the finite dimensional approximation is 

du k = [-u\k\ 2 u k - i ^2 (k ■ u h ) (u\ - ^p' k)] dt + g k d$, kelC N . 

h,ieC]v 

h+l=k 

We set 

Uk = (r k + is k ) J= i, 2 ,3, 

where k • r k = k • s k = and r k , s k , j = 1, 2, 3, are real-valued. Since -u_ k = u^, we are going 
to choose a smaller set of indices k £ JCn in order to take into account that some equations in 
the system are redundant. We set 

K} N = { k £ Z 3 | |kU < N, h > } 

K? N = { k £ Z 3 | |k|oo < N, k 3 = 0, k 2 >0} 

K.% = { k £ Z 3 | |k|oo < N, k 3 = k 2 = 0, h>0} 



and 

in such a way that 



/C = K} N U K? n U JC%, 
)C N = ICU (-£) and K n (-£) = 0. 
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Notice that #(/C) = ^[(2N + 1) 3 — 1], we call such number D. Now, if k e /C, the sum extended 
to all pairs of indices h, 1 such that h + 1 = k can be written in the following way: 

£ = £ + £ + £■ 

h+l=k h+l=k h+l=k h+l=k 

h,ie/Civ h ,ie/c he/c he-C 
ie-/c ie/c 

since if h, 1 ^ /C, k does not belong to /C as well. We denote by \ the sum extended to 
indices in /C. With this position 

(k-0(ui- ^-T^k) = ^*(k.u h )( Ul -^-^k)+ X)^ k ' Uh )(^-7kjr k ) 

h+l=k 1 1 h+l=k 1 1 h-l=k 1 1 

h,lG/Cjv 

(k-u h )(ui- — -k) 

l-h=k ' ' 



so that the equations become 
du k + (u\k\ 2 u k + i ^ (k • u h ) {u\ - ^7j^ k ) + i ^ (k ■ u h ) (^I - TjqT' k ) + 

h+l=k ' ' h-l=k 



+ i (k • w h ) {u\ - -7TTT k )) = 9k 



l-h=k 1 1 

It is convenient to write explicitly the equations relative to the real and imaginary part of u^, 



and 



dry, + (v k 2 r k - 


£*(" 

h+l=k 


■ r h )(si - 




(k 


• s h ) (n - 


|k| 2 j 




+ 


£< k 

h-l=k 


■ r h )(si - 


w k >- 


(k 


■ Sh) (n - 








£*(" 

l-h=k 


■ r h )(si - 


w k )- 


(k 


■ s h ) (n - 




) dt = qldfi 


dsk + (v k 2 s k + 


£>• 

h+l=k 


r h ) (n - 




(k- 


Sh)(si - 


w k > 




+ 


£("■ 

h-l=k 


r h ) (n - 




(k- 


Sh)(si - 


W k) 




+ 


x — 

£(*• 

l-h=k 


r h ) (n - 




(k- 


Sh)(si - 


w k » 


dt = qldft. 
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In view of the above formulas, we set 





F rl = -Hk| 2 ri 


+ £> 

h+l=k 




k • si 

|k| 2 


fci) + (k 




k • n 

|k| 2 


(3.1) 




h-l=k 




k ■ si 

|k| 2 


h) - (k 




k ■ r\ 

|k| 2 






+ 5> 

i — n — k 


n.)W- 


k ■ si 

|k| 2 


fci) - (k 


«h)(r|- 


k • n 

|k| 2 


and 


















F sl = -i/|k|V k 


x — -\* 

- 5>- 

h+l=k 


r h ) (rj - 


k ■ r\ 

|k| 2 


h) -(k- 


«h)(s{- 


k • s\ 

|k| 2 


(3.2) 




- E*(k' 

h-l=k 


r h ) (rj - 


k ■ r\ 

|k| 2 


%i)+(k. 


fh)(si- 


k ■ s\ 

|k| 2 






- £< k - 

l-h=k 


r h ) (rj - 


k ■ r\ 

|k| 2 


**) + (k- 




k ■ s\ 

|k| 2 



h) 
h) 
h). 



4. The Lie Algebra generated by the dynamics 

The state space of the Markov process (r(t),s(t)) which is solution of the equations stated 
above is a linear space U C R 6D , where D = #JC, given by 



(4.1) £/ = 0OR k ©Sk), 



ke/c 

and each element of U is labelled (r, s), with r = (r k , r k , r k ) kg £ and s = (s k , s k , s k ) ke ^, and 

i? k = { (r, s) e R 6D | r k • k = 0, s k = 0, r h = s h = 0, h + k } 
S k = {(r,a) GR 6D |s k -k = 0, r k = 0, r h = s h = 0, h^k}. 

In the same way, we can define the Lie algebra il corresponding to the vector space U, 
(4.2) ii={G\G= J2 ^^ + G sL ^-andk.a k = k.G Sk = 0}. 

k G/C k k 

i=l,2,3 



We define also the subspaces il k = 9^ k © (5 k of il of constant vector fields, where 

= { £ ^in^M and e k = { J3_4^i«k6^} 



i=l,2,3 i=l,2,3 

In this section, we want to find some reasonable conditions on the set M of forced modes 



(such set has been defined in Section 3.1 ) in such a way that the algebra generated by the fields 



(4.3) 



{Fo}uii k keAA 
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where 



kg/C 

1=1,2,3 



and F r i^ and F s i have been defined respectively in (|3.1| ) and Q3.2|) , contains all the constant vector 
fields of il. In particular, it follows that the Hormander condition holds, that is the generated 
Lie algebra, evaluated at each point of U, gives U itself. We start with some computations that 
will be useful in the sequel. 



Lemma 4.1. Let m, n G K and V G il m , W G il n , with 

dr{ 



d 



d 



- + W° , 

dr J m ds J m 



then 



(i) z/k = m + n ; h = n — m and g = m — n ; 
[[^o, V],W] = [(v r ■ k)P k (w s ) + (w s ■ k)P k {v r ) + (v s ■ k)P k {w r ) + {w r ■ k)P k (^' 



_d_ 



+ [ 




k)P k {w s ) 


+ 




k)P k (^ s ) 






■ k)P k ( 


w r ) 






k)p k K 


+ [ 


(v r 


•h)P h «) 


+ 




■h)P h (^ 


— 


(v s 


•h)P h 






(w r 


• h)P h (^ 




(v r 


h)p h K; 


+ 




• h)P h « 


+ 




■h)P h 


« 


+ 


(w s 


• h)P h (i; s 


+ [ 


{v s 


g )p g K) 


+ 


(tu r 


g)^ s ) 


— 


y 






- 








(v r 


■ g)W 


+ 






+ 




S)P 8 ( 




+ 







d 

ds k 

d 
ds h 
d 





ds s 

where P k is the projection of R 3 on the plane orthogonal to the vector k, and in the 
above formula the terms corresponding to indices out of JC are zero; 
(ii) if there is q G Q such that n = qm, then [[P , V], W] = 0, 
(Hi) [[F Q ,V],W] = ±[[F ,V + W],V + W\. 



Proof. We compute the derivatives of the components of F (defined in (|3.1| ) and (|3~2l) ), 



-Hk|%<5 km + kj(s k _ m - s l m _ k + s l m+k ) + k • (s k _ m - s m _ k + s m+k )(^ 



dF r , 

'k 

dr J m 
dF r , 

'k 

8s J m I "I 

^^ s k i I i i i \ l / \ ( Si 2kikj \ 

~= —KjVk-m + r m-k + r m+kJ — k • (Jk-m + r m-k + r m+kj {Oij yT\2~) 



IH2 / 



ki(rl_ m + r. 



'3 v k-m ^ ' m-k ' m+k, 



+ k • (r k _ m + r m _ k - r m+k ) (5,. 



2 /c j 



'j 



ap s i 

A k 



-Z^|k| 5jj5 km + ^j(s k _ m k S m+k) ~l~ ^ ' ( S k-m s m _ k S m _|_ k ) 
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and the second derivatives (we have set A* 7 (k) = Sukj + Sijki — 2 fc ?^j 2 fc ' ), 
d 2 F r , d 2 F r , d 2 F rl 

" - k -0, — — ^ = (£ n ,k-m - ^n,m-k + £ n ,m+k) >&(k) 



dr l n dr J m ds l n dsm ds l n dri 



and 



— (<5 n ,k-m + ^n,m-k + ^n,m+]s.) A 7 (k) , 



. . \~n,is — ill . -11,111 — is. i -11,111-rn./- -o t v— / j , 

or l „orin osLor 



d 2 F si 

~= (^n,k-m — ^n.m-k — ^n,m+k) A\ (k) 



with the agreement that everything concerning indices out of the set /C is zero. 

Let V G il m and W G il n as in the statement of the lemma, then by computing the bracket 
we obtain 

^ Q2 F Q2 F g Q2 F Q2 F g 

We analyse the coefficients of the 5^-components: 

2 F r i d 2 F r i 

v° w\ — -. — \ + v^wf 

<* /JnJ /-Jin/, ■* 



dsindrl drlndsL 



'm 

3 



^(<W-n - 5m,n-k + 5 m ,n+k) -4}j(k)^lt;[ + (<5 n ,k-m ~ <$n,m-k + £n,m+k)^}{ ( k )^J^ 



' 7 <^ 

(<W-n - <W-k + ^m,n+k)[(^ ' k)P k («; r ) i + (w T ■ k) P k (v 

~(^n,k rn ^n,m— k + ^n,m+k) [(^ r ■ k)F k (ty s )j + (w s ■ k)P k (v 



r \), 



where Pk(v)i = v 4 — • k). In a similar way it is possible to treat the coefficients of the 

c^-components, and claim (i) is true. 
If n = q m, it follows that 

v r ■ k = v s • k = w r ■ k = w s ■ k = 

with k = m + n, m — n and n — m, so that using property (i) of this lemma, claim (zz) holds 
true. Finally, if V G il m and G il n , by the Jacobi identity, 

[[F , V + W],V + W] = [[F , V],V} + [[F , V],W} + [[F , W],V] + [[F , W},W} = 2[[F , V],W}. 

□ 

The computations of the above lemma show that the non-linear term mixes and combines the 
components. In some sense, this mechanism can be considered as a geometrical counterpart of 
the cascade of energy. Our aim is to understand for which sets M of forced modes the evaluation 
of the Lie algebra generated by the fields (|4.3|), gives U . We define the set A(Af) C /Cjv of indices 
k G K N such that the constant vector fields corresponding to k (or to — k, depending on k G /C 
or — k G /C) are in the Lie algebra generated by the vector fields (|4.$). Obviously, Af C A(J\f), 
and our aim is to show that A(J\f) = /Cat. 
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Lemma 4.2. Let M be a subset of indices and define the set A{N) as above, 
(i) If me A(Af), then also -m G A(J\f), 

(ii) ifm, n are in A(Af), m + n is in Kn, m and n are linearly independent and |m| ^ |n| ; 
then m + n G A{J\f). 

Proof. The first property follows from the fact that w_ m = u^. In order to show the second 
claim, take m and n in A(J\T) R /C and assume that k = m + n G /C. The claim follows if 
m + n G A(J\f). 
Let 

«=i 111 i=i 111 i=i 11 «=i 11 

with u • m = w ■ n = 0. Then, by property (i) of the previous lemma, 

[[F ,^],W] + [[F ,V],W^] = 2((t;-k)P k H + ( W -k)P k (i;))-A ) 
[[F 0; y],^]-[[F 0) y s ],W s ] = 2(( V -k)P k H + («;.k)P k ^))-A. 

Now, let M, E two vectors in R 3 such that {k, M, E} is a basis of R 3 , M and E span 
{ x G R 3 | x ■ k = } and m, n are in Span[k, M\. Choose 

v = \ x \s. + + v x E, w = A 2 k + fi 2 M + u 2 E, 

then, by the assumptions on m and n, it is always possible to choose the coefficients Ai, /ii, 
v\, A2, (-12, V2 in such a way that (v ■ k)P k (to) + (w ■ k)P k (t;) is any vector in Span[M, E]. 
In other words, il k is contained in the Lie algebra generated by the vector fields ( |4.3| ). In 
the same way, if h = n — m G /C (or if g = m — n G JC), the conclusion follows by taking 
[[F , V r },W s ] - [[F , V s ], W r ] and [[F , V%W r ] + [[F , V%W% □ 

5. Determining sets of indices 
In view of Lemma [4.2|, we call a subset M of /Cjy a determining set of indices for the ultraviolet 



cut-off N, if M generates the cube K N in the sense that A{M) — JCjsr, where A(N) has been 



defined in the previous section. Lemma [4.2| shows us which is the algebraic structure of such 
set. Namely, A(N) is symmetric with respect to the origin and it is close with respect to the 
sum, under some restrictions (m and n have to be linearly independent, with |m| 7^ |n| and 
m+n G /Cjv)- If one neglects such restrictions, Lemma ^72| tells us that a set M is a determining 
set of indices for the cut-off iV if it is an algebraic system of generators for the group (Z 3 , +), 
that is, the smallest subgroup of Z 3 which contains M is the whole Z 3 . 



Since by Lemma [L2]it is obviously true that a determining set of indices, with respect to any 
cut-off, is a system of generators, one can ask if the vice-versa is true, that is if each system of 
generators is a determining set of indices for a suitable cut-off. We give this statement in the 
form of a claim, since in our opinion any proof seems to be full of technicalities which are not 
of great interest in this context. 

Claim 1. If M is an algebraic system of generators for the group (Z 3 , +) and M C JCn, then 
M is a determining set of indices for the ultraviolet cut-off N . 
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For the sake of completeness, we give (see Jacobson j7|, Theorem 3.8 and Theorem 3.9) a 
necessary and sufficient condition for a set of indices Af to be a system of generators of the 
whole group Z 3 . 

Theorem 5.1. A set Af C Z 3 is a system of generators of Z 3 if and only if the g.c.d. of the 

minors of order 3 of the matrix A is equal to 1, where A is the k x 3 matrix whose rows are the 
coordinates of the points of Af and k = #J\[ . 

The intuitive idea which lets us believe that the claim is true is that the restrictions given in 



the statement of property (ii) of Lemma O can be avoided in the following way. 

The restriction about linear independence can be easily avoided by moving aside: for example 
if one wants to sum m with itself, the best way is to obtain 2m as m + n + m — n, where n is 
linear independent with m. 

The restriction about the Euclidean norm (that is, |m| ^ |n|) can be avoided, where possible, 
as in the previous case. Sometimes, as in the case of the proposition below, this is not possible, 
since it may happen that all indices we are allowed to use, have the same Euclidean norm. In 
such a case the solution is to reach the index by different paths, providing with each path a 
component of the Lie algebra we are dealing with, in analogy with Lemma |4.1| . This method is 
probably peculiar of the dimension three and it does not hold in lower dimensions (see E and 
Mattingly ||). 

Indeed these tricks are used in the proof of the following proposition, which states that the 
working example we talked about in Section |2] is a determining set of indices. 

Proposition 5.2. Any set Af C Z 3 containing the three indices (1,0,0), (0,1,0) and (0,0,1) 
is a determining set of indices. 



Proof. A careful analysis of the last part of the proof of Lemma |4.2| shows that, if |m| = |n|, then 
the Lie brackets [[F , V], W], with V G ilm and W G il n , span the two-dimensional subspace of 
iWn given by 

d o 
XE-- + fJ,E-- , 

where E is the index orthogonal (in R 3 ) to m and n. Hence, if we sum (1,0,0) and (0, 1,0), 
we obtain the corresponding two dimensional subspace of 11(1,1,0). Again a direct computation 
shows that such a smaller subspace is indeed sufficient, since if we combine it with il(o,o,i) we 
obtain the full 11(1,1,1). Now we just subtract (0, 0, 1) from (1, 1, 1) to obtain the full il(i,i,o) and, 
in the same way, we can obtain all the indices of norm 2. With this set of indices is now easy 
to obtain, by means of Lemma [4.2| and of the tricks explained above, all the indices in JCn, 



whatever is N. □ 

As a consequence of the above proposition is that if Af is a determining set of indices for a 
cut-off N, then it is a determining set of indices for any other cut-off threshold larger than N. 



6. The control problem 

The section is devoted to the proof of the controllability properties of the finite dimensional 
approximations of Navier-Stokes equations. The first part contains some generalities on poly- 
nomial control systems. The approach and the results are taken from Jurdjevic and Kupka 
0. In the second part we adapt the proof of a theorem (again of Jurdjevic and Kupka ||) 
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to our case. The original theorem applied to polynomials of odd degree. Polynomials of even 
degree behave in a different way, mostly because of the obstructions of the positive terms. Our 
case has no obstructions, essentially because of property (ii) of Lemma |4.1| , and the system is 
controllable. 

6.1. Generalities on polynomial control systems. We consider a system of the form 

m 

X = p(x) + 

8=1 

where x G R n , &i, b 2 , . . . , b n are fixed vectors in R n and P is a polynomial mapping, that is 
P = [Pi, . . . , P n ) and each Pi is a polynomial in the variables (xi, . . . , x n ). Let Yi, . . . , Y n be 
the constant vector fields assuming respectively value b\ , . . . , b n and let F be the vector fields 
having the components of P as its components, and define 

F = { F + ^2 u iYi | • • • , u m ) G R m } 

i=l m 

We define, for each xq G R™ and t > 0, the set y4jr(x ,t) of states reachable, with a suitable 
control it = (u\, . . . , u n ), from the initial state xq in a time smaller than t. We define the set 
A^(x , t) of states reachable exactly at time t. 

Two families of vector fields T\, Ti are said to be equivalent if for all x G R n and t > 0, 

AfiOM) = Af 2 (x,t). 

If JF is equivalent to T\ and to JF 2 , then it is equivalent to T\ U JF 2 . It makes sense then to 
define the saturate of J 7 , denoted by Sat (J 7 ), which is the union of all families of vector fields 
equivalent to T . Moreover we will call Lie(JF) the Lie algebra generated by T . Finally, the Lie 
saturate of T is defined as LS(JF) = Sat(jF) n Lie(jF). In order to obtain controllability, the 
Lie saturate should be as large as possible, as stated by the following theorem. 

Theorem 6.1. Let T be any family of smooth vector fields and assume that LS(JF) contains n 
vectors Vx, . . . , V n such that the vector space spanned by them is in LS(JF) and for each x G R n 
the vectors Vi(x), . . . , V n (x) span R n . Then A^(x, t) = R n for each x G R ra and t > 0. 

We adapt the conclusions of the theorem to the system that will be studied in the following 
section. 

Corollary 6.2. Let J 7 be any family of smooth vector fields and assume that the constant vector 
fields o/LSfT) span R n . Then A*/(x,t) = R n for each x G R n and t > 0. 

Proof. From the previous theorem, Aj?(x, t) = R n . Moreover, by Theorem 13, Chapter 3 of 
Jurdjevic (see also the remarks after Theorem 11, Chapter 5 of ||) it follows that also 
A*p(x,t) = R n . □ 

In the following, we will need the following two lemmata. The first lemma permits the 
enlargement of a family of vector fields by means of diffeomorphisms. A diffeomorphism (f> : 
R n — > R n is a normaliser of a family T if for all x G R n and t > 0, 

4>(Ar(4>-\x),t)) c A r (x,t), 

we will denote by Norm(jF) the set of all smooth normaliser of T . 
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Lemma 6.3. The family T is equivalent to \J ( p &iorra ^{4>*{y) \ V G J 7 }, where 0* is the differ- 
ential of (j). 

The second lemma gives the geometrical structure of the Lie saturate of a family of vector 
fields. 

Lemma 6.4. // T is any family of smooth vector fields, then T is equivalent to the closed 
convex cone generated by {XV | < A < 1, V <E J 7 }, where the closure is in the C°° topology on 
compact sets of~R n . 

6.2. Control of the finite dimensional approximations of Navier-Stokes. We are able 
now to prove the controllability property of our equations. We aim to prove that the control 
problem 



s k -F Sk (r,s) =gjX> 



(6.1) 

where F rk and F Sk are defined in ( p.l|) and (|3~2|), and the 3x3 matrices are defined in (|3.1D, is 
controllable, in the sense that for each initial state (r/, si) G U, for each final state (rp,SF) G U 
and for each time T > there is a family of controls (t£, vQ^^f, where Af is the set of indices 
corresponding to the non-zero q^, such that the solution corresponding to that control starts 
at t = in (rj, si) and arrives in (r>, sf) at time t = T. 

Theorem 6.5. Assume that the set Af of non-zero components of the control is a determining 
set of indices, as defined in Section [|. Then system (|6.1|) is controllable in the sense given 
above. 

Proof. First we show that ilk C LS(JF) for k G Af. Let A G R and k G Af and take Yk G ilk, 
since 

AYk = lim — (F + nAF k )) 

n^oo n 

it follows by Lemma |6.4j that Yk £ LS(JF). 



Now we aim to show the following claim: if m, n G JC^ are linear independent indices with 
|m| 7^ |n| and m + n G /Cat, and if il m © il n C LS(JF), then also il m+n C LS(JF). If the claim 
is true, it follows that each ilk is contained in LS(JF) and, since by the assumptions Af is a 
determining set of indices, Corollary |672| applies and the proof is ended. 



The proof of the claim now follows. From Lemma |4.2| we know that il m+n is spanned by 



[[Fq, V), W], where V G il m and W G il n . From property (in) of Lemma [O], we have that 



[[Fo, V},W} = i[[F , V + W], V + W] and - [[F , V], W] = ^[[F , V - W],V - W]. 



Since by Lemma [BT^ LS(JF) is convex, in order to prove the claim it is sufficient to show that 
X[[F , V),V] is in LS^) for each A > and V G il m © il n . 

So, let a G R and V G ilm©il n , then <f){x) = e aV (x) is in Norm(jF) (see the proof of Theorem 
2 of Jurdjevic and Kupka 0), so that, by Lemma |6~3| , (e ay ) ! „(Fo) G LS(JF). Now, since the 
coefficients of F are polynomials of degree 2 and V is a constant vector field, it follows that 

(e aV UF ) = I + a[V, F ] + -[V, [V, F \], 
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and so, for each A > 0, 

\[V,[V,F )) = hm ^(e aV )*(F ) E LS(F), 

a^oo a 

since LS(JF) is closed. The theorem is proved. □ 

7. The exponential convergence 



In this last section we prove Theorem |2.2| as a consequence of a general result by Meyn and 
Tweedie [0] (see Theorem 6.1). Before giving the statement of such theorem, we need to state 
some definitions. They will be given in a simplified form, adapted to our case, while the general 
statements can be found in the papers by Meyn and Tweedie ||11|| , [12]. 

A nonempty subset C of the state space U is a petite set for a Markov process with transition 
probabilities -Pt(v) if there are a non-trivial measure ip and a probability distribution a on 
(0, oo) such that 



J Pt(x, •) a(dt) > tp for all x G C. 



A function V : U — > R + is a Lyapunov function for the process if V(x) — >• oo as \x\ — > oo 
and there are real constants c > and d such that 

CV(x) < -cV(x) + d 

where £ is the generator of the diffusion. 
The kinetic energy 



i 2 , j 2 X 

>k 

will play the role of the Lyapunov function in our case, as stated by the following lemma. 



Lemma 7.1. For each (r,s) 6 £7, 

3 

^^(r^ 4 + 4£ 4 )=0, 



kec 



i=l 



where the polynomial E^ and E s i^ are respectively the homogeneous part of degree 2 of the 
polynomials F r i and F s i , defined in (|3.1[ ) and (13.21), and 

CV(r,s) < -2isV(r,s) + o- 2 , 
where £ is the generator defined in ( |2.2| ) and a 2 zs the variance of the Brownian motion B t . 



Proof. The first property is a far consequence of the famous property of the non-linear part of 
Navier-Stokes equations, namely J v ■ (v ■ V)n = 0. Indeed 

(7.1) ^2( r k ■ E rk + s k ■ E Sk ) = 5^ Im [( k ' M h)(wi • Uk)] = Uk ' Eu * 

keic keic h,ieic N keK N 

h+i=k 

where u k = r^ + isk and E Uk is the non-linear part in equation fl2.1|), namely 



h,l6/Cjv 
h+l=k 
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Finally, the proof that the last sum in (|7. 1| ) is equal to is just a matter of swapping the two 
indices k and 1. 

The second property is then an easy consequence of the previous one: 

CV = ^(-2z/|k| 2 (r£ + 4) + r k • + s k • £, k ) + a 2 < -2vV + a 2 

□ 

Now we are able to prove Theorem 



Proof of Theorem From Theorem [2.1| we know that the Markov process (r(t), s(t)) is strong 
Feller and irreducible. Using Theorem 3.3 and Theorem 4.1 of Meyn and Tweedie ||11|| , it follows 
that all compact sets of the state space U are petite sets. Moreover the previous lemma tells 
us that the kinetic energy V is a Lyapunov function. By means of Theorem 6.1 of Meyn 
and Tweedie [0, we conclude that there are positive constants C and p such that for each 
(r ,s ) e U, 



a 2 



\\Pt((r , s ), •) - Try, < C(l + V(r , s ) + — Je"*. 
with f = l + V. □ 
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